Experimental improvement of chaotic synchronization due to 
multiplicative time— correlated Gaussian noise 
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The effect of time-correlated zero-mean Gaussian noise on chaotic synchronization is analyzed 
experimentally in small-size arrays of Chua's circuits. Depending on the correlation time, an im- 
provement of the synchronization is found for different values of the noise amplitude and coupling 
diffusion between circuits. 



Recently there has been considerable interest in 
stochastic resonance, i.e., the enhanced response of a sys- 
tem to an external signal induced by noise [Wiesenfeld 
& Moss, 1995; Gammaitoni, et al., 1998; Luchinsky, et 
al., 1998], a phenomenon in which noise has a creative 
role. Moreover, the influence of noise has also been stud- 
ied within the context of arrays of cells [Lindner, et al., 
1996; Braiman, et al., 1995a;1995b; Gailey, et al., 1997; 
Shuai, et al., 1998]. Here, the coupling strength and the 
noise intensity play an important role for array enhanced 
stochastic resonance [Lindner, et al., 1995]. Applications 
of noise to biological systems or in engineering problems 
could be of special relevance. On the other hand, the 
behavior of uncoupled chaotic systems under the influ- 
ence of external noise has been the subject of recent work 
[Maritan & Banavar, 1994; Pikovsky, et al., 1994; Herzel 
& Freund, 1995; Malescio, 1996; Gade & Basu, 1996; 
Longa, et al., 1996; Shinbrot, et al., 1993; Sanchez, et 
al., 1997;1999]. The main idea behind these papers is 
that uncoupled chaotic systems cannot be synchronized 
by means of an identical noise signal (Gaussian noise of 
zero mean), except for a noise with some non-zero bias. 

In this Letter, the role of a time correlated Gaussian 
noise on diffusively chaotic coupled cells is analyzed. The 
dynamical noise used in this Letter is a Gaussian noise 
of zero mean of the Ornstein-Uhlenbeck type [Sancho, et 
al., 1982], characterized by a correlation function 
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where r is the correlation time and a is the noise am- 
plitude. In the limit r — > the white-noise limit is 
recovered. We emphasize here the case of global noise, 
where the noise is identical at each site, as opposed to 
the case of incoherent or local noise, where the noise is 
uncorrelated from site to site. 

Experiments have been performed with N resistively 
coupled Chua's circuits [N — 3) [Madan, 1993; Chua, 
1995] in the chaotic regime, accordingly to the design in- 
troduced by [Sanchez, et al., 1997;1999]. Each circuit, 
j = 1,2,3, in the array is defined by the following evolu- 
tion equations, 
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where V\, Vi, and the voltages across C\ and C-2 and 
the current through L, respectively, are the three vari- 
ables that describe the dynamical system, resulting from 
a straightforward application of Kirchhoff 's law. The pa- 
rameters have the following meaning: C\ and Ci are the 
two capacitances, L the inductance, R the resistance that 
couples the two capacitors and ro the inner resistance of 
the inductor. Circuits were connected through capacitor 
C\ by resistances R c , leading to a diffusion term in the 
potential differences [Chua, 1995], with a coupling coef- 
ficient D cx 1/R C . Circuits at the boundaries are only 
connected with one neighbor. 

The three-segment piecewise-linear characteristic of 
the nonlinear resistor (Chua's diode) is defined by, 

h(Vi) = G b V 1 + ^(G a -G b )[\V 1 +B p \-\V 1 - B p \] (3) 

where G a and G b are the slopes of the inner and outer 
regions of h(V\), respectively, and B p = 1 V defines the 
location of the breaking points of the three-slope nonlin- 
ear characteristic h(V\). 

An experimental setup of three identical Chua's cir- 
cuits driven by noise has been built. Their components 
are defined by 

(C 1 ,C 2 ,L,r ,R) = (lQnF,100nF,10mH,20fi,l.lkQ). 
The slopes of the nonlinear characteristic h(V), Eq. (^), 
are defined by G a = -8/7000 and G b = -5/7000. Fig- 
ure ^ shows the schematic diagram of the experimental 
setup. The circuits were sampled with a digital oscillo- 
scope (Hewlett-Packard 54825A) with a maximum sam- 
pling rate of 4 x 10 9 samples per second, 1.5 GHz band- 
with, and a record length of 32000 points, connected to 
a PC for data processing. 

The external noise has been introduced multiplica- 
tively using a recently introduced circuit [Sanchez, et al., 
1997;1999] that enables to drive the nonlinear element by 
using the voltage from an external source. The nonlin- 
ear element is driven, in general, by the voltage coming 
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from an external source, not necessarily the voltage com- 
ing from capacitor C±, as may happen in the case of a 
standard Chua's circuit. Thus, it is a voltage controlled 
current source (VCCS) with a characteristic defined by 
Eq. (||). This yields the following evolution equation for 
the voltage across the capacitor Ci, 

CxV 1 = ^^-h{V 1 +m) (4) 

where it is easy to see that the noise term yields a mul- 
tiplicative contribution. 

The time correlated noise, (Eq. (|l|)), is obtained 
electronically by passing the output voltage of a white 
Gaussian noise generator, £ w (t), through a single-pole 
active filter with a time constant r = RbCf,, before being 
applied to the circuit as it is shown in Fig. [l| [McClintock 
& Moss, 1989; Luchinsky, et al., 1998]. The external 
white noise has been generated by using a function gen- 
erator (Hewlett-Packard 33 120 A). Their characteristics, 
gaussian distribution and zero mean in the absence of an 
offset, have been adequately checked^. 

Figures ^(a-d) show the physical properties of an ex- 
perimental time-correlated noise before being added to 
the voltage V\ in the Chua's circuit. The power spec- 
trum of the noise (b) cannot be considered to be flat 
within the frequency range of interest, r" 1 w 200 Hz. For 
small noise intensities, the obtained probability distribu- 
tion does not fit perfectly to the Gaussian distribution as 
a consequence of the experimental noise in the sampling 
process (Fig. ||(c)). The correlation function was also cal- 
culated experimentally (d) and compared with the theo- 
retical function given by Eq. ([!]) . Our realization of noise 
decays exponentially within a time scale almost equal to 
the expected theoretical value. 

In order to characterize the degree of synchronization 
between cells of the array, we introduce the following 
time-averaged quantity, 

t=l \ 3 =2 J 

This function is positive defined and it is equal to zero 
when all the cells in the array are globally synchronized. 
Since K can serve as a measure of the array complexity, 
in this context it can be related to the Kolmogorov-Sinai 
entropy [Benettin, et al., 1976; Klimontovich, 1996]. 

The main effect of a colored Gaussian noise on an ar- 
ray of diffusively coupled chaotic cells is to improve the 



*As the bandwith (« 10 MHz) of the white noise from the 
HP generator is much higher than the characteristic frequen- 
cies of the Chua's circuit, then for our purposes, we can con- 
sider this ideal "white" noise. 



synchronization between units. Figure |3J shows the evo- 
lution of if as a function of the correlation time r for 
different values of the coupling resistance R c . The mean 
value of K increases with R c as expected, so a scaling 
factor was introduced in order to compare the different 
observed behaviors of K. In general, independently of 
the specific value of R c as r increases, first K decreases 
almost exponentially, reaches a minimum, and then rises 
smoothly until a constant saturation value is attained for 
t > 1, The minimum of K {K m i n ) corresponds to an 
optimum choice of r (r m j n ) to obtain the best synchro- 
nization. 

For intermediate values of r, the time correlated Gaus- 
sian noise periodically modulates V\ when driving the 
nonlinear element. A resonance effect between the 
Chua's time scale and the noise correlation time, r, 
should be expected, since the power spectrum of the 
noise cannot be considered to be flat within the frequency 
range of interest, r _1 . This resonance effect could explain 
the improvement of chaotic synchronization observed in 
Fig. ||for K = K min [Lorenzo & Perez-Munuzuri , 1999]. 
Here, the double-scroll attractor becomes periodically 
asymmetric with increasing noise amplitude, as well as 
blurred (Figs. ||(a-b)), due to the slow dynamics of the 
noise. 

The value of r corresponding to K m i n was found to 
increase with the coupling resistance as it is shown in 
Fig. [|. As t — > T m i n , a stronger interaction between 
the two characteristic time scales of both cell and noise 
should be expected, then improving the synchronization 
between circuits. For r = T m j„, in terms of frequency 
locking, the dynamics of the cell could be simplified to 
that of an oscillator forced periodically with a frequency 
equal to . Then, for a chain of linearly coupled oscil- 
lators, its dynamics can be described in terms of a plane 
wave solution with a wave velocity proportional to \[T), 
The wave dispersion relation is given by u> cx yD/X, with 
u) the wave frequency and A the wave length. For small 
size arrays, it can be considered that A is fixed by the 
boundary conditions. In this case, the wave frequency 
increases with the coupling strength. Thus, in order to 
obtain locking between the internal oscillation frequency 
and the external periodic forcing, as the coupling resis- 
tance increases, the external forcing period (in our case 
this is related to the value of T m i n ) should also increase. 
Obviously, the explanation above is a simplification of the 
problem, since the chaotic dynamics cannot be mapped 
in a simple way to that of an oscillator. Nevertheless, 
our aim is to stress the similarity between the classical 
frequency locking problem that occurs in a chain of os- 
cillators forced periodically and the behavior of K for 
t « T m i n . Here, the locking does not occur for a single 
value of the frequency, but for a range of frequencies that 
gives rise to a wide behavior of K as a function of r near 
the onset of resonance. 

On the other hand, the two limits, r — > and r — > oo, 
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in Fig. |3| deserve further comments. When r — > 0, the 
white Gaussian noise limit is recovered and circuits in 
the array do not become synchronized to each other in- 
dependently of the variance of the noise [Sanchez, et al., 
1997;1999]. In this case, the structure of the unperturbed 
double-scroll attractor gets smeared out with the increas- 
ing noise amplitude, corresponding to a decreasing signal- 
to-noise ratio, while no evidence of synchronization be- 
havior is observed. Similarly, when r — > oo the term 
in Eq. (^) behaves as a constant value added to the volt- 
age V\ . Noise affects the double-scroll dynamics that be- 
comes asymmetric, while no synchronization is observed 
between cells within the array. In fact, for high enough 
noise amplitude, the main effect will be a biased signal 
that will induce a regularization in the system. This ef- 
fect is analogous to that of some chaos suppression meth- 
ods that achieve this result through perturbations in the 
system variables [Matfas & Giiemez, 1994; 1996]. 

Figure || shows the dependence of K with the noise 
amplitude for a small value of r near the white limit 
case. As expected, increasing the noise strength leads to 
a worse synchronization between circuits (K increases). 

The influence of the number of circuits in the array, 
N, on chaotic synchronization by time correlated noise 
is also being studied [Lorenzo & Perez-Munuzuri, 1999]. 
Preliminar results show that the larger the array, the 
larger the needed fluctuations to improve chaotic syn- 
chronization; i.e. the variance of the noise must increase 
as N increases. 

We have observed a non-monotonic dependence of the 
degree of synchronization measured in terms of K as a 
function of the noise correlation time r and the coupling 
between cells in a one-dimensional array. In fact, for val- 
ues of t of the order of the time scale of the chaotic attra- 
tor a stochastic resonance effect is found that could ex- 
plain the observed minimum value of K . In other words, 
the effect of color noise to improve the chaotic synchro- 
nization between cells is more robust than that of white 
noise for a constant noise amplitude. 
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FIG. 1. Schematic diagram of the experimental setup used 
to introduce noise in a multiplicative way in an array of dif- 
fusively coupled Chua's circuits. Noise is added to voltage 
Vi and is used to drive the nonlinear element of the circuit 
(see Eq. (§)). The noise is buffered from all VCCS's ensur- 
ing thereby no interaction between the circuits except that 
due to the coupling resistances R c . The white Gaussian noise 
generated by the function generator is transformed by the 
high-cut filter in a Gaussian noise of zero mean of the Orn- 
stein-Uhlenbeck type (see Eq. (|l])) with a correlation time 
r = RbCt, which is then added to the signal from capacitor 
C\ (also buffered). The low-frequency cutoff of the filter, de- 
termined by R a C a , is fixed at « 1 Hz, and the high-frequency 
cutoff is adjustable by tuning Rb (Ct = ljuF). The output 
noise finally passes through a variable-gain operational am- 
plifier (not shown) before being applied to the circuits. 



FIG. 2. Characterization of a time-correlated Gaussian 
noise obtained after a white zero-mean Gaussian noise is 
passed through a single-pole filter, (a) Temporal evolution 
of the color noise with amplitude 500 mV (peak-to-peak) and 
correlation time r = 5 ms. (b) Power spectrum of the noise 
signal, (c) Probability function distribution. The line shows 
a fitting of the experimental results to a Gaussian curve, (d) 
Correlation function (continuous line). The dashed line cor- 
responds to the theoretical curve given by Eq. (hi). 



FIG. 3. Log-linear plot of K as a function of the time cor- 
relation r for three different values of the coupling resistance 
R c . Since the length of the time series is determined by the 
record length of the oscilloscope and the time scale of the 
Chua's circuits, in order to improve the statistic, 50 realiza- 
tions of the noise were carried out for each value of r. In do- 
ing so, we substitute the limiting value of T — » oo, in Eq. (|^), 
for T > r, in such a way that the most probable value of 
K becomes identical to the ensemble average of K when the 
number of experiments (realizations of noise) is large. The 
experimental data are shown as symbols, while lines repre- 
sent an interpolation to the previous ones: crosses (x) and a 
dashed line for R c = 6.8 Q, rhombi (O) and a dot-dashed line 
for R c = 14 Q, and circles (o) and a solid line for R c = 27 fl. 
The obtained values of K were scaled between and 1 for a 
better representation of the phenomenon. The minimum and 
maximum values of K for each coupling resistance are: 0.11 
and 0.13 for R c = 6.8 Q, 0.34 and 0.37 for for R c = 14 Q, and 
1.39 and 1.51 for R c = 27 Q. Noise amplitude (peak-to-peak): 
250 mV. 

FIG. 4. Effect of a time-correlated noise on the dou- 
ble-scroll (chaotic) attractor (a) of the Chua's circuit for in- 
termediate values of r. As r increases, the attractor becomes 
smeared out as well it becomes "periodically" asymmetric, 
finally loosing the double-scroll appearence (b). The image 
presented in panel (b) was acquired after the oscilloscope was 
stopped at a given instant of time. Otherwise, other possible 
asymmetric shapes of the attractor could have been obtained. 
Noise amplitude (peak-to-peak): (a) 10 mV and (b) 400 mV. 
Correlation time; r = 50 ms. R c — 6.8 Q. 

FIG. 5. Dependence of the time correlation value, corre- 
sponding to the minimum of K, with the coupling resistance, 
R c . Noise amplitude (peak-to-peak): 250 mV. 

FIG. 6. Dependence of K with the noise amplitude 
(peak-to-peak) of a time correlated noise. Note the linear 
dependence observed for the range of used noise amplitudes. 
R c = WD., t = 0.01 s. 
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